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Dijet production is studied in collisions of quasi-real photons radiated by the LEP 
beams at e+e“ centre-of-mass energies \/Sgg = 161 and 172 GeV. Jets are recon¬ 
structed using a cone jet finding algorithm in the range | < 2 and Ey* > 3 GeV. 
The angular distributions of direct and double-resolved processes and the inclusive 
two-jet cross-section are measured. They are compared to next-to-leading order 
perturbative QGD calculations and the prediction of the leading order Monte Garlo 
generators PYTHIA and PHOJET. 


1 Introduction 

The production of dijet events in the collision of two quasi-real photons is 
used to study the structure of the photon and different QCD predictions. The 
measurement of inclusive jet cross-sections and the comparison with next-to- 
leading order (NLO) QCD calculations and different photon structure functions 
can constrain the relatively unknown gluonic content of the photon. 

2 Event selection and jet finding 

To select a data sample of two-photon events the following set of cuts was 
applied. The sum of all energy deposits in the electromagnetic calorimeter 
(ECAL) and the hadronic calorimeter (HCAL) has to be less than 45 GeV. 
The visible invariant hadronic mass, Wecal, measured in the ECAL has to be 
greater than 3 GeV. The missing transverse energy of the event measured in 
the ECAL and the forward calorimeters has to be less than 5 GeV. At least 5 
tracks must have been found in the tracking chambers. Events with detected 
scattered electrons (single-tagged or double-tagged events) are excluded from 
the analysis. 

The results of the cone jet finding algorithm depend on the minimal trans¬ 
verse energy E™'" and the cone size R = with the pseudo¬ 

rapidity r] = — lntan(0/2) and the azimuthal angle cj). Here the values were 
chosen to be i? = 1 and E™'" = 3 GeV. The jet pseudorapidity in the labora¬ 
tory frame is required to be within < 2. 
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After applying all cuts and requiring at least two jets 2681 events remain 
in the data corresponding to an integrated luminosity of 20 pb“^. For the 
simulation of two-photon interactions the Monte Carlo generators PYTHIAEI 
and PHOJET □ have been used. The mean of the selected Monte Carlo 
events is 0.06 GeV^. About 1.2 % of the events in the data sample are expected 
to be e+e" annihilation events with hadronic final states and 0.2 % electron- 
photon events. 


3 Angular distributions in direct and resolved events 


A pair of variables, and , can be defined which is related to the fraction 
of the photon energy participating in the hard scattering: 


(E+Pz) 


I jets=l,2 
= 


(E+Pz) 


hadrons 


and 


jets—1,2 

E (E-Pz)’ 

hadrons 


( 1 ) 


where Pz is the momentum component along the 2 axis of the detector and E 
is the energy of the jets or hadrons. If an event contains more than two jets, 
the two jets with the highest Ei^^ values are taken. Samples with large direct 
and double-resolved contributions can be separated by requiring both xji and 
x~ to be larger than 0.8 (denoted as x^ > Q.8) or both values to be smaller 
than 0.8 (denoted as x^ < 0.8), respectively^. 

The transverse energy flow around the jet axis for double-resolved events is 
expected to show additional activity outside the jet due to the photon remnants 
compared to the direct events. Figure ^ therefore shows the transverse energy 
flow around the jets with respect to the jet direction for data samples with 
different x^ cuts. The pseudorapidity difference is defined by: 

Af]' = k{f] — r?-’®*). 

The factor k is chosen event-by-event to be /c = -1-1 for events with > x~ 
and k = —1 for events with xjl < x~. Due to the definition of Aij' there is 
always more of the remnant found at Ay' < 0 and the enhancement due to the 
additional transverse energy flow at negative and positive Ay' is asymmetric. 
The jets in the data sample with x^ > 0.8 are more collimated and there is 
almost no activity outside the jet whereas the transverse energy flow of two-jet 
events with x^ < 0.8 shows additional activity outside the jets. 

In the centre-of-mass system of the interacting partons or bare photons 
the parton scattering angle 9* is defined as the angle between the jet axis of 
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Figure 1: Uncorrected energy flow transverse to the beam direction measured relative to the 
direction of each jet in two-jet events. Jets from data samples with a large contribution of 
(a) double-resolved and (b) direct events according to their a;ijr and Xj values are shown. 
The energy flow is integrated over |A<^| < 7r/2. Statistical errors only are shown. 


the jets originating from the outgoing partons and the axis of the incoming 
partons or bare photons. The angular distribution of the jets can be used 
to check the separation between the different event classes using a variable 
which is not directly correlated to the definition of and x~. In the dijet 
centre-of-mass frame one expects different angular distributions for direct and 
double-resolved events. An estimator of the angle between the jets and the 
parton-parton axis in the dijet centre-of-mass frame can be formed from their 
pseudorapidities. The variable cos0* is calculated as 


cos 9* = tanh 



Since the ordering of the jets is arbitrary, only | cos 9* \ can be measured. The 
matrix elen^nts of elastic parton-parton scattering processes have been calcu¬ 
lated in LOq. The cut on * restricts the accessible range of values of | cos 9* \. 
Requiring the invariant mass of the dijet system to be larger than 12 GeV en¬ 
sures that values of | cos0*| < 0.85 are not biased by the cut. The boost 
of the two-jet system in the z direction is defined hy fj = ( 77 !®*^ -|- 
The detector resolution on | cos0*| deteriorates significantly for events with |? 7 | 
larger than 1. These events were therefore rejected by requiring \fj\ < 1. 

Figure || shows the |cos0* [-distribution of events with x^ > 0.8 and of 
events with x^ < 0.8. The dependence on the Monte Carlo models used 
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Figure 2: Angular distribution of events with large direct and large double-resolved con¬ 
tributions according to the reparation with Xj and . The data are compared to QCD 
matrix element calculationsu. The data are normalised to have an average value of f in the 
first three bins and the curves are normalised to be 1 at cos(O). 


is taken into account by adding the difference between the results obtained 
with PYTHIA, which are taken to be the central values, and PHOJET to the 
systematic error. The error bars show the statistical and the systematic errors 
added in quadrature. The events with > 0.8 show a small rise with | cos6>*|, 
whereas the events with < 0.8 show a much stronger rise in | cos0*|. The 
data points of the events with x^ < 0.8 are compared with the prediction of a 
QCD matrix element calculation of the interaction of quarks or gluons in the 
photon. The matrix elements with a relevant contribution to the cross-section 
where antiquarks are involved instead of quarks show a similar behaviour as 
the examples shown. The QCD matrix element calculations agree well with 
the data points of the data samples with large direct and large double-resolved 
contribution. 


4 Inclusive two-jet cross-sections 

The inclusive two-jet cross-section is measured using a cone jet finding algo¬ 
rithm with a cone size R = 1. The data were corrected for the selection cuts, 
the resolution effects of the detector and the background from non-signal pro¬ 
cesses. In Fig. the inclusive two-jet cross-section is shown as a function of 
Ay*. The bin sizes, which are indicated by the vertical lines at the top of the 
figure, approximately reflect the experimental resolution. 
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Figure 3: The inclusive two-jet cross-section as a function oLF?^* for jets with |? 7 j®*| < 2 
compared to the NLO calculation by Kleinwort and Kramer 0. The direct, single-resolved 
and double-resolved cross-sections and the sum (continuous line) are shown separately. 

The -By* distribution is compared to an NLO perturbative QCD calcula¬ 
tion of the inclusive two-jet cross-section by Kleinwort and KrarnerS who use 
the NLO GRV parametrisation of the photon structure function □. Their cal¬ 
culation was repeated for the kinematic conditions of this analysis. The direct, 
single- and double-resolved parts and their sum are shown separately. The 
resolved cross-sections dominate in the region By* ^ 8 GeV, whereas, at high 
By* the direct cross-section is largest. The data points are in good agreement 
with the calculation except in the first bin where the NLO calculation is not 
reliable due to IR singularities. The uncertainties due to the modelling of the 
hadronisation process are expected to contribute mainly at low By* values. 

The inclusive two-jet cross-section, which is dominated by the low By* 
events, depends on the parton density function of the photon which mainly 
differs in the assumptions on the gluon part of the photon. In 40 to 60 % 
of all processes at least one gluon from the photon interacts depending on 
the parton density function used. This leads to different predictions of the 
inclusive two-jet cross-section. The inclusive two-jet cross-section as a func¬ 
tion of |r?^®*| is shown in Fig. |[ The inclusive two-jet cross-section predicted 
by the two Monte Garlo models differ significantly even if the same photon 
structure function is used. This model dependence currently reduces the sen¬ 
sitivity to the parametrisatioii-|Of the photon structure function. The GRV-LO 
and SaS-lD parametrisationsHQ describe the data equally well, but the LAGl 
parametrisationS overestimates the inclusive two-jet cross-section significantly. 









FREIBURG-EHEP-97-05 


6 


^ 800 
S700 
^600 
I 500 
400 
300 
200 
100 
0 


- 1 -r 

OPAL preliminary 

jetl iet2 ^ 

> 3 GeV 


- 1 - 

PYTHIA LACl 
PHOJET LACl 
PYTHIA GRV 
PHOJET GRV 
PYTHIA SaS-ID 


0 0.2 0.4 0.6 0.8 


1.2 1.4 1.6 1.8 


h' 


jet. 


Figure 4: The inclusive two-jet cross-section as a function of | for jets with > 3 GeV. 


5 Conclusions 

The distribution of the parton scattering angle 0* of data samples with large 
direct and double-resolved contributions separated experimentally using the 
variables and x~ have been compared to the relevant QCD matrix element 
calculations. The inclusive two-jet cross-sections were measured as a function 
of Ei^^ and The Ei^^ dependent two-jet cross-section is in good agree¬ 

ment with an NLO QCD calculation. The GRV-LO and SaS-lD parametri- 
sations describe the inclusive two-jet cross-section equally well. The LACl 
parametrisation overestimates the inclusive two-jet cross-section significantly. 
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